Nuclear multifragmentation time-scale and fluctuations of largest fragment size 
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We show that the biggest fragment charge distribution in central collisions of ^■^^Xe+"''*Sn leading 
to multifragmentation is an admixture of two asymptotic distributions observed for the lowest and 
highest bombarding energies. The evolution of the relative weights of the two components with 
bombarding energy is shown to be analogous to that observed as a function of time for the largest 
cluster produced in irreversible fragmentation for a finite system. This analogy allows us to speak 
of pseudo-critical nuclear multifragmentation around Ebeam ~ 30 MeV/A. We infer that the size 
distribution of the largest fragment in nuclear multifragmentation is also characteristic of the time 
scale of the process, which is largely determined by the onset of radial expansion in this energy 
range. 

PACS numbers: 25.70.Pq 



Phase transitions occupy a central place in many fields 
of physics such as condensed matter, collisions of atomic 
aggregates, and most recently, in formation of quark- 
gluon plasma in ultra-relativistic hadron collisions. They 
allow to investigate the equation of state and phase di- 
agram of the system under study. The case of nuclear 
multifragmentation, as observed in intermediate energy 
heavy- ion collisions provides a unique opportunity 
to study not only thermodynamical properties of nuclear 
matter, but also phase transitions in finite systems. 

The universal character of order parameter fluctua- 
tions in finite systems [l, Q , provides a good framework 
in which to address such questions. Using such a model- 
independent analysis, it was shown that the size (atomic 
number) of the largest fragment produced in multifrag- 
mentation events, Zmax, behaves like an order parameter, 
i.e. the scaling properties of its fluctuations change with 
increasing energy In this sense nuclear multifrag- 
mentation belongs to a large class of phase transitions 
or critical phenomena for which the size of the largest 
cluster is the order parameter, for example site- or bond- 
percolation on lattices, reversible or irreversible aggrega- 
tion . Using the latter class of model we will show that 
we can extend this analogy, by studying the evolution of 
the order parameter distribution in the critical region. 

In it was shown that a distinct asymptotic form of 
the largest fragment distribution in multifragmentation 



can be associated with each scaling regime : a quasi- 
Gaussian symmetric Zmax distribution at low bombard- 
ing energies, and an asymmetric Gumbel distribution Q 
in the high-energy disordered regime. In this letter we 
will study in more detail the transition from one regime 
to the other, using new data on Z^ax distributions for 
i23Xe-|-"°*Sn central collisions measured with INDRA 
at bombarding energies intermediate between the two 
asymptotic regimes. We will show that, in fact, the 
Zmax distributions evolve continuously from one asymp- 
totic distribution to the other with increasing bombard- 
ing energy, and, by analogy with the evolution of the 
order parameter distribution in the irreversible aggrega- 
tion model, deduce a link between Z^ax fluctuations and 
the time-scale of nuclear multifragmentation. 

The irreversible aggregation model describes an out-of- 
equilibrium clusterization process, whose order parame- 
ter is the average size of the largest cluster at any given 
time, Sniax- The evolution of the droplet size population 
is given by the Smoluchowski equations Q , 

A'yQCj ^^K,,c^c, (1) 

i-\-j—s i—1 

where Cs is the concentration of droplets of size s, Kap = 
K[)af3 is the coalescence probability per unit of time be- 
tween two droplets of respective size a and /?, Kq is a 
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Figure 1: Standard deviation of the largest cluster size, 
o"(smax), as a function of the time t, for irreversible aggre- 
gation with A*' = 216; t*, t~ and t"*" are the pseudo-critical 
time and the limits of the critical domain (see text). 



constant and N is the size of the system. This model 
exhibits a second-order phase transition after a critical 
time, tc 0. 

Precise limits of the critical domain for finite systems, 
where the correlation length becomes comparable to the 
size of the system Q , are not easily obtained. An estima- 
tion using the width of the largest cluster size distribu- 
tion, fT(sinax)i h^s bccu proposcd in In the infinite 
system, this critical quantity must diverge at the criti- 
cal time. As shown in Fig. 1 [2 for a finite system with 
N = 216 particles, cr(sniax) reaches a sharp maximum at 
the pseudo-critical time, t* = 1.28tc- We define the lower 
and upper limits of the critical domain, and i+ , as the 
times at which (T(smax) is above 10% of its maximum: 
t- = OMc and t+ = 2Mc (see Fig. [J). 

The order parameter distribution for irreversible aggre- 
gation is not known exactly, but it was proposed that the 
asymptotic Smax distribution for t <JC is that of Gumbel 
[10| . The average size of the largest cluster evolves with 
time as more and more coalescence of smaller clusters 
takes place. At long times, the order parameter is then 
essentially of additive nature. From the central limit the- 
orem, this results in an asymptotic Gaussian distribution 
of Smax for i » t * . In the critical domain, for t~ < t < , 
fluctuations are so large that similarly prepared systems 
can exhibit one or the other behavior. This suggests that 
the effective Smax distribution for an ensemble of sys- 
tems, at any given intermediate time, is an admixture of 
a Gaussian distribution and a Gumbel distribution. 

We performed numerical simulations of the Smolu- 
chowski equations, Eq.([T]), for finite systems with ^ 
50 — 200 particles, comparable with nuclear system sizes. 
The obtained s^ax distributions were then fitted with 
the following admixture of Gaussian and Gumbel distri- 



butions: 
where, 
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and X =Sniax- The five parameters (77, a, a,n and 
arc left to vary freely, apart from the constraint that r] 
should be between zero and one. Fits were performed 
using a binned maximum-likelihood estimation method 



11| . In order to quantify the evolution of the relative 



weight of each distribution, we introduce the ratio: 



R = -J — —j— = 2?? - 



1 



(3) 



where Ica and Igu are the integral of the components de- 
fined in Eq.dH). This ratio equals -1-1 for a pure Gaussian 
distribution and —1 for a pure Gumbel distribution. 

The results of the fitting procedure applied to the dis- 
tribution of the largest cluster size between t~ and 
are presented in Fig. [5] At short times, t~ < t < t^, the 
distribution is almost purely "Gumbellian", with a small 
Gaussian component (Fig. [2ja)), and the ratio i? « — 1 
(Fig. [2l^e)). As time increases, the Gaussian compo- 
nent increases, such that between the infinite-system crit- 
ical time tc (Fig. [2l^b)) and the pseudo-critical time 
(Fig. djc)), both components have comparable weights, 
and the ratio R passes through zero (Fig. [UJe)). The 
very large fluctuations of Smax observed at these times 
in FigH] are due to the superposition of the two distri- 
butions, whose combined width is far greater than the 
intrinsic variance of each component. Finally, for the 
longest time-scales, the Gumbel component disappears 
(Fig. WLd)) and R tends towards -1-1 (Fig. ^e)). 

The analysis of Smoluchowski calculations shows that 
the decomposition of the Smax distributions into an ad- 
mixture of two asymptotic distributions seems to be rel- 
evant. Similar results were obtained for all simulated 
system sizes, N: as N increases the transition becomes 
more abrupt and the value R = occurs closer and closer 
to tc- It should be noted that we have obtained similar 
results with a three-dimensional percolation model [l^ . 
which describes the at-equilibrium formation of clusters 
on a geometric lattice as a function of bond probability, 
p, with a second-order transition at the critical proba- 
bility, pc. Such models have often been used to study 
criticality in nuclear multifragmentation [13l |14| | . In all 
cases the value R = occurs in the critical region of the 
underlying phase transition, between the infinite-system 
critical point and the A'^-dependent maximum of Smax 
fluctuations. We will now apply the same analysis to 
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Figure 2: (color online) Analysis of Smoluchowski calculations 
for TV = 216 particles: (a-d) distribution of the size of the 
largest cluster P(smax) = (l/A^)dA''/dsmax for different times 
in the critical domain; (Black solid curve) best fit to the data 
using Eq.(2); (Red dashed curve) Gumbel component; (Blue 
dotted curve) Gaussian component; (e) evolution of the ratio 
K between the two components (Eq.(3)) in the critical domain 
(see Fig. [T]), as a function of the time t in unit of tc (see text); 
vertical lines are the same as in Fig. [1] 



the largest fragment distributions measured for nuclear 
multifragmentation in central heavy-ion collisions. 

Collisions of ^^^Xe-|-"°*Sn were measured using the IN- 
DRA 47r array at the GANIL accelerator facility. This 
charged product detector, composed of 336 detection cells 
arranged according to 17 rings centered on the beam axis, 
covers 90% of the solid arigle. More experimental de- 
tails can be found in 0, [la | . Data were taken during 
two separate campaigns. Beam energies of 27, 29, and 
35 MeV/A have been performed specifically in order to 
probe in more detail the energy range around the change 
of scaling regime observed in [J] : those data are presented 
here for the first time. 

As in previous works [1, @| j we study multifragmenta- 
tion events occurring in central collisions, requiring the 
geometrical overlap between projectile and target to be 
as close as possible to total. In addition, the correct mea- 
surement of the largest fragment distribution requires the 
detection of nearly all fragments event by event. For this 
purpose we will use two global variables, Eti2 and Ztot- 



Figure 3: (color online) (Points) Log-log correlation between 
the first two cumulant moments ((Zmax)^and cr^) of the order 
parameter's distribution; (Lines) linear fits performed in the 
range 50-32 MeV/A (resp. 32-25 MeV/A) which correspond 
to a slope A ~ 1 (resp. A ^ 0.5) (see text). 



Eti2 is the total transverse energy of light charged par- 
ticles {Z ~ 1,2), which increases with collision centrality 
[lij , while Ztot is the sum of the atomic numbers of all 
detected charged products in each event. In sepa- 
rate energy-dependent cuts were made in the Eti2 and 
Ztot distributions. In this new analysis, in order to en- 
sure equivalent selections for all beam energies we retain 
events which maximize the quantity (Eti2xZtot)- For 
each beam energy, a cut was defined corresponding to 
the last percentile of the (Eti2xZtot) distribution mea- 
sured with the online trigger condition (charged product 
multiphcity M > 4). 

Let us first examine the scaling properties of Zmax fluc- 
tuations including the three new data points at 27, 29 
and 35 MeV/A. Figure [3] shows the correlations between 
the first two cumulant moments of the Zmax distribution 
((Zmax)^and cr^), for each beam energy. As in Fig. 1(c) 
of the data fall on two branches, ^ (Zmai}^^, 
with different values for the scaling parameter, A [^Q: 
A - 1 above 32 MeV/A, and A - 1/2 below 32 MeV/A. 
The new data points present a consistent behavior which 
follows the systematic scaling trend. It should also be 
noted that the scaling properties of Zmax distributions 
in central collisions are quite independent of the detailed 
selection procedure, as different methods were used to 
select events in 0, Q . 

The results of the fits performed on the Z^ax distri- 
bution using Eq.([2|) are shown in FigH) It can be seen 
that, for all analyzed energies, a good fit to the exper- 
imental distribution is achieved: reduced values all 
lie between approximately 4 and 10, and can further be 
reduced if Z^ax distributions are smoothed to remove 
odd-even staggering of the yields [17| : in that case all 
values are close to 2, except for data at 32 MeV/A for 
which a significantly larger value (x^ ~ 4) is obtained. 



At the highest energy considered here (50 McV/A), the 
result of the fit is an almost pure Gumbel distribution 
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Figure 4: (color online) Analysis of ^■^^Xe+"'''Sn reactions: 
(a-d) distribution of the charge of the largest detected frag- 
ment P(Zniax) = {l/N)dN/dZnia.x for different beam energies, 
(Black solid curve) best fit to the data using Eq.(2), (Red 
dashed curve) Gumbel component, (Blue dotted curve) Gaus- 
sian component; and (e) evolution of the ratio R between the 
two components (Eq.(3)), as a function of the beam energy. 



(Fig. mjd)). For the other bombarding energies, both 
contributions arc present, and the relative importance of 
the Gaussian component increases with decreasing en- 
ergy. The quantity R, defined in Eq.(l3]), is evaluated for 
each bombarding energy, and its evolution is presented 
in Fig. m^e). The vertical bars show the uncertainty 
coming from the fitting procedure: (t_r = 2(7,;, calculated 
from Eq. ([5]) . The value i? = is reached between 29 and 
32 MeV/A, in the same bombarding energy range as the 
change of A-scaling (Fig. [3|). 

A strong similarity can be seen between the results 
of the analysis for irreversible aggregation (Fig. [2]) and 
central collisions of ^^^Xe+"'^*Sn (Fig. |4]). In the aggre- 
gation model the order parameter distribution depends 
on the length of time during which clusters can form. 
In central heavy-ion collisions, a determining factor for 
the time-scale of fragment formation is the radial ex- 
pansion of the multifragmenting system, which increases 
with bombarding energy [ij. Fragment sizes can evolve 
as long as exchanges of nucleons can take place between 
them, i.e. until freeze-out at which the short-range nu- 
clear interaction between fragments ceases and they be- 



gin to move apart. It has been shown that for central 
^^^Xe-|-"°*Sn reactions the onset of significant radial ex- 
pansion occurs for beam energies above 25 MeV/A[l5|. 
The similarity between Figs. [5] and SI can therefore be 
understood in terms of fragment size distributions be- 
ing determined on shorter and shorter time-scales due to 
increasingly rapid expansion. 

We can use this interpretation to understand the sys- 
tematic evolution of A-scaling with beam energy and 
system mass presented in Q: for ^^Ni-|-^®Ni collisions, 
measured from 32 MeV/A to 90 MeV/A, a change of 
A-scaling and of the form of the Z^ax distribution were 
observed, as for -'^^^Xe-|-"''*Sn but at a higher bombarding 
energy of 52 MeV/A; for the lighter system ^'^Ar+KCl, 
studied from 32 MeV/A to 74 MeV/A, only the A = 
1/2 regime was observed, with quasi-Gaussian Zmax dis- 
tributions; on the other hand, for the much heavier 
^^^Au-f^^^Au system, at bombarding energies between 
40 MeV/A and 80 MeV/A, only the A = 1 regime oc- 
curs, with Gumbel Zmax distributions. 

Radial expansion in central heavy-ion collisions oc- 
curs after significant compression of the incoming nu- 
clear fluid, and as such depends not only on static nu- 
clear matter properties such as incompressibility, but also 
on transport properties such as the degree of stopping 
achieved in the collision. The latter increases with the 
mass of the colliding nuclei, as shown in [l^]- Thus 
for light systems, such as ^'^Ar-l-KCl or ^^Ni-|-^^Ni, the 
bombarding energy required to achieve sufficient initial 
compression for there to be significant radial expansion 
is higher than for the heavier systems like ^^^Xe-|-""*Sn 
and ^^^Au-|-^^^Au. This explains why the A-scaling tran- 
sition occurs at higher energy for ^^Ni-|-^^Ni than for 
i29Xe-h"°*Sn . For the very Hght ^^^Ar-^KCl system the 
threshold is higher still than for ^^Ni-|-^®Ni, outside the 
range of measured bombarding energies. On the other 
hand, for -'^^^Au-l-^'^^^Au both the greater initial compres- 
sion and the far larger Coulomb contribution may come 
into play in order to reduce the fragment formation time- 
scale even at the lowest bombarding energy. 

We have shown that, for finite systems, the largest 
cluster size distribution in critical aggregation models 
is an admixture of the two asymptotic distributions ob- 
served far below and above the critical region. This result 
holds true for both equilibrium (percolation) and out-of- 
equilibrium (irreversible aggregation) models. A similar 
decomposition has been shown for the experimentally- 
observed charge distribution of the largest fragment per 
event produced in nuclear multifragmentation, indicating 
that the critical domain lies around Ebeam ~ 30MeV/A 
for the ^^^Xe-l-^^'Sn system. By analogy with the irre- 
versible aggregation model, where the form of the order 
parameter distribution depends on the time-scale of the 
process, we interpret such criticality along with the cor- 
responding change of A-scaling as the onset of an "explo- 
sive" multifragmentation regime in which initially com- 
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pressed heated nuclear liquid clusterizes in the presence 
of significant radial expansion [l9| . The mass-dependence 
of the energy at which the onset occurs for different (sym- 
metric) systems is related to nuclear stopping and hence 
to transport properties of hot nuclear matter. Such an 
overall picture is both consistent with, and provides a link 
between, recent results on the role of radial expansion in 
nuclear multifragmentation 



15, m. 



and the systematic 
study of nuclear stopping around the Fermi energy [l^l ■ 
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